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We investigate the effect of next nearest-neighbour hopping on the zero temperature Drude weight 
or superfluidity in mesoscopic one-dimensional systems of (a) single particle with quasidisorder 
(Aubry-Andre model) and (b) hardcore bosons with nearest-neighbour interaction. We show that 
there is an interaction-induced enhancement of the Drude weight when the next nearest-neighbour 
hopping is frustrated for the many-body system. The observed non-monotonic behaviour of the 
Drude weight occurs because the repulsive interactions first suppress the frustration in the system, 
leading to a rise in the Drude weight, whereas at much larger interaction strengths the charge density 
wave fluctuations set in the insulator and the Drude weight drops. The present work reveals a novel 
scenario in which a persistent flow enhancement is plausible in the presence of kinetic frustration 
and repulsive interactions. 
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I. INTRODUCTION 

An important feature distinguishing quantum me¬ 
chanical systems from classical systems is the Aharonov- 
Bohm effect [T|, wherein a charged particle picks up a 
phase while moving through a field-free region which en¬ 
closes a finite magnetic flux. Such a flux is proportional 
to the vector potential along the path of the particle and 
perturbs the energy levels of the quantum mechanical 
system [H . A notable example of this effect occurs in low 
temperature mesoscopic rings in the presence of a static 
magnetic field, where the electronic wave function can 
retain the Aharonov-Bohm phase acquired throughout 
the ring and coherently extend over the whole system, 
thereby generating a persistent current @1 Hi- This is 
true even if the material is not superconducting, as had 
been predicted by theoretical studies and corroborated 
by experiments on copper Q and gold @ rings, although 
these persistent currents are orders of magnitude smaller 
than those in superconductors Q . Remarkably persistent 
currents can also be generated in systems of neutral par¬ 
ticles, where an artificial gauge field or a rotation of 
the system causes the particles to acquire an Aharonov- 
Bohm phase. Indeed, a persistent flow has been achieved 
in an ultracold setup of Bose-Einstein condensed (BEC) 
atoms trapped in a toroidal geometry, by transferrin g or - 
bital angular momentum to the atoms of the system [lOj . 

A primary question in past studies, prompted mainly 
by the disagreement between experimental findings in 
small metallic rings and single-particle calculations, has 
been to discern the influence of interactions and disorder 
on the persistent currents in one- and two-dimensional 
model systems, see e.g. Refs. [iMl and references 
therein. It is generally understood that repulsive inter¬ 
actions can counteract the localising effects of disorder, 
thereby increasing the persistent current in the system 


[IMl. However, such an interaction-induced enhance¬ 
ment of the current is by no means predicted in all model 
systems (T 3 |. or it is expected only to a weak extent 
[la [l^ . With the observation of a BEC superflow and 
the ability to artificially gauge model Hamiltonians, ul¬ 
tracold atomic setups seem ideal to investigate the effects 
of interactions and disorder on the persistent current in 
various physical models. 

In this work we investigate the zero-temperature 
Drude weight or conductivity stiffness (which determines 
the persistent response to an external field, and which 
is precisely defined below) in two prototypical systems 
that have been realised in cold atom experiments: single¬ 
particle model with incommensurate potential (Aubry- 
Andre model) and hardcore bosons model [isl - l^ at 
and away from half-filling. In particular, we ask what 
the effect of a longer-range hopping of the particles will 
be on the properties of each of these models defined on 
ring geometries. Our motivation for this specific question 
is two-fold: (i) the capability to g enerate artificial gauge 
fields [^, persistent flows |lClj |. and tunable nearest- 
and next nearest-neighbour hopping amplitudes 2 l| in 
ultracold systems raises the question of how longer range 
interactions and hopping processes of the gauged parti¬ 
cles will affect the experimentally observable results; (ii) 
ladder materials, to which our long-range hopping one¬ 
dimensional systems are equivalent, realise rich physics 
through a variety of spin models[ 22 ||. Moreover, their 
bosonic analogues have been studied in the literature and 
novel phases such as a supersolid have been recently pre¬ 
dicted m, [25. 

We now briefly summarise previous related work and 
results. Although the problem of a longer-range hopping 
has been addressed before for the aforementioned models 
[H, [ 25 l - [^ , these works focused largely on critical, gap 
and localisation properties. We will uncover some novel 



2 


aspects missed in these earlier studies, namely that the 
interplay between kinetic frustration and non-local repul¬ 
sive interactions produces a non-monotonic behaviour of 
the Drude weight as the interaction strength is increased. 
Our results will bear similarities to earlier reports on 
frustration-induced metallicity undertaken in the 

context of analysing the phases of the transition metal ox¬ 
ide PrBa 2 Cu 408 and of two-dimensional organic conduc¬ 
tors. In these studies (i) the Ising gap was seen to vanish 
for a range of frustrated interaction strengths in a spin- 
system on a ladder geometry [2^ and such a behaviour 
was interpreted as the onset of a conducting phase due 
to frustration. However, no analysis of the conductivity 
within the metallic phase was attempted; (ii) the Drude 
weight was seen to be enhanced by frustrated interactions 
in a quarter-filled two-dimensional electronic system . 
Our work extends and complements these studies by con¬ 
sidering a kinetically-frustrated bosonic system where we 
find the Drude weight or superfluid density (and hence 
the persistent flow) to be enhanced by repulsive interac¬ 
tions, even in the thermodynamic limit; that is, in the 
presence of frustration interactions can indeed enhance 
conductivity. We also explicitly demonstrate how such an 
effect is absent both for a single-particle one-dimensional 
system which undergoes a genuine conductor-insulator 
transition, and for the many-bod y sy stem in the absence 
of frustration. A recent study of interacting one¬ 
dimensional bosons in the continuum with a localised and 
moying barrier in a ring reached a qualitatiyely similar 
conclusion: that there is an optimal on-site interaction 
for which the persistent flow is maximal. 

Drude weight and superfluid density 

The persistent flow in a normal metal, superconduc¬ 
tor or superfluid may be induced by threading a flux 
through a toroidal or ring system. Here we consider a 
ring of L sites threaded by a flux cf. This flux giyes rise 
to a static yector potential along the path of the parti¬ 
cles, which modifies the hopping amplitudes by a Peierls 
phase factor [3l| {t —>■ and perturbs the energy 

levels of the quantum mechanical system Q . This sensi¬ 
tivity of the ground-state energy E to the external flux 
produces a thermodynamic persistent flow I oc —dE/dcj) 
around the ring [3, . 

The response of the ground-state energy to an in¬ 
finitesimal flux is measured by the Drude weight ttD, 
which quantifies the strength of the zero-frequency peak 
in the real part of the Kubo conductiyity. In one¬ 
dimensional systems it is giyen by Issl - ls^ 

L{E{(j)) - E{(j)o)) ^ 1 

{(j) — (foY 2 \m) ’ 

( 1 ) 

where 0o is the flux at which the ground-state energy is 
minimum and {p/m) is the effective ratio of the den¬ 
sity of mobile carriers to mass. Thus, for an infinitesimal 
flux, the Drude weight provides a way of quantifying the 
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FIG. 1: (Colour online) Renormalised Drude weight D (see 
text) in extended Aubry-Andre model as a function of the 
quasi-commensurate potential strength V. 


persistent flow in the ring m- Equivalently, the Drude 
weight can be understood as a measure of the sensitiv¬ 
ity of the ground state to a small twist cj) in the bound¬ 
ary conditions. Consequently, a vanishing (finite) Drude 
weight signals an insulating (conducting) phase at zero 
temperature [Ij. Finally, we point out that the superfluid 
density is defined very similarly as the Drude weight in 
eq. (P). Specifically, it corresponds to calculating first 
the ground-state energy in the thermodynamic limit and 
then its second derivative with respect to the flux [s^. 
Nonetheless, in (quasi) one-dimensional systems D and 
the superfluid density are equal due to the finite num¬ 
ber of energy level crossings in the thermodynamic limit 
[a Hi. Therefore, our results for the Drude weight may 
also be interpreted as a measure of the superfluid density 
in the system. 

In what follows we compute, using exact diagonal- 
isation, the ground-state energy of the two aforemen¬ 
tioned models (both with and without a flux) defined 
on chains of length L with periodic boundary conditions. 
The largest matrix we diagonalise has a linear dimension 
of about 10 millon. 


II. SINGLE PARTICLE SYSTEM 

In this section we extend the results presented in Ref. 
[a |- regards the superfluid density or Drude weight 
of the Aubry-Andre model [33]) to the case with next 
nearest-neighbour (NNN) hopping. The model, defined 
on a ring of L sites, reads 


L d^E{(f) 

2 d(Y 
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L L 

H = -ti + h.c) - t 2 '^{b]bj ^2 + ^i-c) 

j 3 

L 

j 

where j runs over the lattice sites, b (6'!') is the standard 
annihilation (creation) operator and n is the correspond¬ 
ing particle number operator. The hopping amplitudes 
to the nearest and next nearest neighbours are set re¬ 
spectively by ti and t 2 , and ej = V cos(2t: jg) with V the 
potential strength and g = (\/5 + l)/2 the golden ratio. 
We note that the sign of <i in eq. ([U does not have any 
effect on the physics of the model. On the other hand, 
a negative value of t 2 does not allow the energy asso¬ 
ciated with nearest and next-nearest neighbour hopping 
processes to be simultaneously minimized i.e. it yields 
kinetic frustration. 

For a single particle with no incommensurate poten¬ 
tial (V = 0), eq. (P) may be readily evaluated to give 
D = (ti + At 2 )/L = t^s/L = p/ 2771, where p and m 
are the system’s density and the particle’s bare mass. 
In order to facilitate comparison of results between dif¬ 
ferent ^2 values, we define a renormalised Drude weight 
D ■= LD/{tes) which yields £) = 1 for a single particle 
in the absence of disorder for any value of t 2 ■ 

We consider a chain length of Fibonacci number 
L = 610 in order to impose periodic boundary condi¬ 
tions [ 2 ^, and compute the Drude weight for the single¬ 
particle Aubry-Andre model with NNN hopping ampli¬ 
tudes t 2 = 0, ±0.2 ti. Our results are shown in Fig.Pand 
bear out the following expectations: (i) there is a gen¬ 
uine conductor-insulator transition as a function of the 
potential strength V in the ground state, as seen by the 
sharp drop in the Drude weights. Indeed, the transition 
point V = 2 ti is known exactly for the nearest-neighbour 
Aubry-Andre model; (ii) a negative value of the hopping 
t 2 yields a reduced tes and works as a localising mecha¬ 
nism. Consequently, a weaker quasi-disorder is sufficient 
to fully localise the wave function. In contrast, the crit¬ 
ical disorder strength increases for a positive ^ 2 - These 
results are in quantitative agreement with inverse partici¬ 
pation ratio (IPR) and Shannon entropy calculations . 
However, we note that through our Drude weight calcu¬ 
lations the conductor-insulator transition in the ground 
state may be pinpointed more accurately than through 
IPR calculations. The primary result of this section — 
that the localising effect of the incommensurate poten¬ 
tial always produces a monotonic decay of the Drude 
weight, although much slower for the unfrustrated sys¬ 
tem — will be contrasted with the delocalising effect of 
the non-local repulsive interaction (and its corresponding 
effect on the persistent flow) while competing with ki¬ 
netic frustration: in particular, the competition between 
interaction and frustration will be shown to produce a 
non-monotonic decrease of the Drude weight before the 


conductor-insulator transition sets in. 

III. HARDCORE BOSON SYSTEMS 

Next we consider a system of hardcore bosons de¬ 
scribed by the Hamiltonian 

^hch — H ± D ^ , (3) 

i 

where V now sets the strength of nearest-neighbour inter¬ 
action, ej = 0 and the particle operators satisfy the usual 
hardcore boson commutation relations. Such a model has 
been extensively investigated recently [2i,|2| in the con¬ 
text of analysing the competition between quantum fluc¬ 
tuations and kinetic frustration induced by the ^2 hop¬ 
ping term, and whether such a scenario can lead to exotic 
phases. We will work along a contour of the phase dia¬ 
gram presented by Mishra et al. (Fig. 2 in Ref. SI) 
and uncover some unique properties of the system in the 
conducting phase. 

Given that the second derivative in eq. © is eval¬ 
uated at the location /jq where the ground-state energy 
E{(j)) is minimum, the dependence of the many-body en¬ 
ergy levels on the flux deserves some explanation. In the 
absence of next nearest-neighbour hopping, the hardcore 
boson model is exactly solvable and the Drude weight 
may be written in closed form [33| . The model may also 
be Jordan-Wigner transformed to a system of spinless 
electrons and most results from the bosonic model trans¬ 
fer to the fermionic one. However, when a flux (j) threads 
the ring of electrons the energy is minimum either at 
(/ = 0 or TT, depending on whether the number N of elec¬ 
trons on the ring is odd or even, respectively. This shift 
in energy is due to the statistical phase that arises when 
an electron goes around the ring, thereby changing its 
place with the other N —I electrons. Such a parity effect 
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FIG. 2: Energy levels as a function of flux in a 16-site 
ring with 2 hardcore bosons and system parameters t 2 /ti = 
0.2,y/ti = 1. 
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does not arise in the corresponding hardcore bosonic sys¬ 
tem whose many-body wave function is symmetric under 
particle interchange. Thus, the energy minimum is at 
(j)Q = Q and this is the point about which eq. o is eval¬ 
uated. The location of the energy minimum when NNN 
hopping is included has also been found to be at (po = 0, 
provided the system is not in the incommensurately or¬ 
dered phase obtained in Ref. [2^. In Fig. [5] we show 
a representative plot of the flux dependence of the first 
three energy levels, for a system of two hardcore bosons 
over a 16-site ring, with t 2 = 0.2 ti and V = ti. The en¬ 
ergy crossing at (/) = tt results in a jump discontinuity of 
the persistent flow, and the energy levels — as well as all 
physical properties — will clearly be periodic functions of 
(j) with period 27r. 
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Conductor-insulator transition 


A transition from a conducting to an insulating 
(charge density wave) phase is expected as the nearest- 
neighbour repulsive interaction strength V is increased, 
and this is true whether t 2 is present or not [H, [s^. 
In the latter case, the Berezinskii-Kosterlitz-Thouless 
(BKT) transition point is known to be exactly atV = 2 ti 
for a half-filled system. On the other hand, in the pres¬ 
ence of NNN hopping, the BKT transition points have 
been computed numerically by Mishra et al. |23l | from the 
behaviour of the single particle excitation gap, using the 
density matrix renormalisation group technique. Here we 
focus on the conducting phase of the system at half-filling 
and consider the representative frustrated value of 0 = 
—0.2ti for which the phase transition was predicted 23| 
to occur at 14 ~ 1.4 ti. Our main results are summarised 
in Fig.|3]for a r ang e of system lengths. As expected from 
previous works [J[3 IH, HI] , we find the Drude weight 
to scale polynomially D{L) = D + a/L + h/Lp' in the con¬ 
ducting phase and exponentially D{L) = C'exp(—L/^) 
in the insulating phase, where ^ is the localisation length 
which diverges at the transition, as shown in the bottom 
panel of Fig. [3] 

We point out that we do not, nonetheless, employ 
the scaling laws to identify the transition point. Instead, 
we assume the transition takes place at 14 = 1-4 ti [3| 
and perform a polynomial extrapolation to the thermo¬ 
dynamic limit of the computed Drude weights within the 
conducting phase. Such an extrapolation, shown as the 
full line in the top panel of Fig. [3l works well both for 
t 2 = — 0.2ti and for the exactly solvable case of ^2 = 0 
(inset of figure, showing a g:ree ment of the extrapolation 
with Bethe Ansatz results [3^). For a finite t 2 a similar 
sharp jump in the Drude weight should occur as L ^ oo, 
which expectation is borne out from the plots. 

We now expatiate on the singular interplay of in¬ 
teractions and frustration, and its effect on the Drude 
weight. The most important feature of the results 
shown in Fig. |3] is the rise of the Drude weight as the 
nearest-neighbour interaction is increased, which trans¬ 
lates as an enhancement of the persistent flow in the 



FIG. 3: (Colour online) Drude weight of half-filled bosonic 
rings for various chain lengths L. Top panel: Non¬ 
monotonicity of the Drude weight with increasing interac¬ 
tion for frustrated nearest-neighbour hopping t 2 = —0.2ti, 
and polynomial extrapolation (full lines) from data for L = 
10,12, ... 26 in the conducting phase. Inset shows results from 
a similar extrapolation for t 2 = 0 compared to Bethe Ansatz 
(BA) results. Bottom panel: Scaling of the Drude weight for 
t 2 = 0 (triangles) and t 2 = —0.2ti (circles) as a function of 
inverse system length for interaction strengths V = 1 (open 
symbols) and F = 4 (closed symbols). Full lines are expo¬ 
nential fits D{L) = C'exp(—L/^) and dashed lines are poly¬ 
nomial fits D{L) = D + ajL -I- hjl? to the thermodynamic 
limit L —^ oo. The hts for F = 4 give a localisation length 
^ « 15,11 for t 2 = 0, —0.2ti respectively. As F decreases the 
fits give increasing estimates for the localisation lengths 4 


ring. As pointed out in the introduction, an interaction- 
induced enhancement of persistent flows and conducting 
behaviour has been reported in earlier theoretical studies 
of strongly correlated electronic [lTI - [3 . [^ . hard-J3 and 
soft-core bosonic ^ systems, in different [lll - ll3j and 
similar [3 - 1^ scenarios. Here we uncover a mechanism 
which produces an enhancement of the Drude weight in 
ladder geometries as a consequence of the competition 
between non-local repulsive interactions and kinetic frus- 
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tration within the conducting phase, even after a finite- 
size scaling to the thermodynamic limit is performed. 
Indeed, kinetic frustration favours a situation in which 
the hardcore bosons gain kinetic energy by hopping back 
and forth between a pair of nearest-nei ghb our sites and 
decoupling from the rest of the system [^. The repul¬ 
sive nearest-neighbour interaction V, on the other hand, 
favours a charge density wave like ordering where the par¬ 
ticles tend to avoid sitting next to each other. Thus, ini¬ 
tially increasing V causes a spread of the particles and a 
concomitant increase in charge density wave fluctuations, 
which in turn reduces the effective frustration, thereby 
increasing the Drude weight. However, upon further in¬ 
crease of V the insulating behaviour dominates and the 
Drude weight drops accordingly. And crucially, as seen 
in Fig. 131 such a scenario of interaction-induced enhance¬ 
ment of D is valid for finite rings as well as in the ther¬ 
modynamic limit. 

This result is in fact more generic and we have 
checked its validity away from half-filling. The main dif¬ 
ference in that case being a much more gradual decay 
of the Drude weight after it reaches its maximum, which 
is presumably due to the absence of an insulating phase 
away from half-filling. We have also simulated a system 
with unfrustrated hopping (t 2 > 0) and found a mono¬ 
tonic decay of D with increasing interaction, thereby con¬ 
firming the crucial role played by the frustration: in the 
absence of kinetic frustration only one localising mecha¬ 
nism is present. 

IV. CONCLUSIONS 

We have studied the interplay of kinetic frustration 
with quasi-disorder and interactions on the Drude weight 
of single-particle and many-body (quasi) one-dimensional 
systems, respectively. For the single-particle Aubry- 
Andre model, we found that the Drude weight drops 
monotonically with an increase in the strength of the 
incommensurate potential and that, in agreement with 
earlier work [^ . kinetic frustration allows for a quicker 


onset of the delocalisation-localisation transition. 

For the many-body system of nearest-neighbour 
hardcore interacting bosons, we found that a frustrat¬ 
ing nearest neighbour hopping causes an interaction- 
induced enhancement of the Drude weight in the con¬ 
ducting phase. This effect arises because of the degra¬ 
dation of the frustration by the nearest-neighbour inter¬ 
action, causing an increase of the Drude weight; as the 
non-local repulsive interactions are further increased the 
insulating behaviour sets in and the Drude weight drops. 
Our study on kinetically frustrated hardcore bosons is 
complementary to an earlier study of a ID spin system 

where frustrated interactions was seen to destroy the 
single-particle gap, and a later related work [2^ which 
showed the enhancement of the Drude weight by frus¬ 
trated interactions in quarter-filled two-dimensional elec¬ 
tronic systems. 

Our work therefore provides evidence for a unique 
situation where the conductivity is enhanced in simple 
kinetically frustrated bosonic systems, and opens up the 
possibility of future cold-atom experiments and theoret¬ 
ical studies wherein the competition between interaction 
and frustration could produce novel effects in the Drude 
conductivity and persistent flows, even in the simplest of 
(quasi) one-dimensional systems. In fact, studies and 
experimental realisations of interacting and frustrated 
systems in low dimensions are many and varied in the 
present day. Moreover, there still remains much to be 
gleaned about the additional effect of disorder on the 
conductivity of these kinetically frustrated systems, and 
on the effects of such a competition in the predicted su¬ 
persolid phase [ 2 ^. 
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